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ABSTRACT 

We present the first blind intcrferomctric detection and imaging of a millisecond radio transient 
with an observation of transient pulsar J0628+0909. We developed a special observing mode of the 
Karl G. Jansky Very Large Array (VLA) to produce correlated data products (i.e., visibilities and 
images) on a time scale of 10 ms. Correlated data effectively produce thousands of beams on the sky 
that can localize sources anywhere over a wide field of view. We used this new observing mode to find 
and image pulses from the rotating radio transient (RRAT) J0628+0909, improving its localization by 
two orders of magnitude. Since the location of the RRAT was only approximately known when first 
observed, we searched for transients using a wide-field detection algorithm based on the bispectrum, an 
interferometric closure quantity. Over 16 minutes of observing, this algorithm detected one transient 
offset roughly 1' from its nominal location; this allowed us to image the RRAT to localize it with an 
accuracy of l'/6. With a priori knowledge of the RRAT location, a traditional beamforming search of 
the same data found two, lower significance pulses. The refined RRAT position excludes all potential 
multiwavelength counterparts, limiting its optical luminosity to < 1.1 x 10 31 erg s _1 and excluding 
its association with a young, luminous neutron star. 

Subject headings: pulsars: individual (J0628+0909), stars: neutron, techniques: intcrferomctric 



1. INTRODUCTION 

The fast (< 1 s) radio transient radio sky is pop- 
ulated by a wide range of source classes and phys- 
ical processes. Neutron star pulses can be as brief 
as a nanosecond, making them the brightest ob- 
jects in the universe and driving our understanding 
of coherent emission (Hankins ct al. 2003). Obser- 
vations at fast time scales are sensitive to pulsars 
and magnetars (Camilo et al. 2000), cxtragalactic tran- 
sients (McLaughlin & Cordes 2003; Lorimcr et al. 2007), 
stars (Osten & Bastian 2008), and planets (Farrell et al. 
1999). Propagation effects, like dispersion and scatter- 
ing, are also detectable at GHz frequencies on millisecond 
time scales and can be used to probe the interstellar and 
intergalactic media (Cordes & Lazio 2002). 

A new class of fast transient is the rotating radio 
transient (RRAT; McLaughlin et al. 2006). Like pul- 
sars, RRATs are rotating neutron stars that emit ra- 
dio pulses. Unlike pulsars, RRATs pulse rarely (roughly 
once out of every thousand rotations) , so they are only 
detectable by their individual pulses. That observa- 
tional distinction suggests that RRATs arc physically 
distinct from ordinary pulsars (although that is not yet 
clear; Keane et al. 2011). Identifying optical and X- 
ray counterparts to radio-detected RRATs would allow 
us to study the parent neutron star and define their 
relation to ordinary pulsars (Kondratiev et al. 2009). 
Of the few dozen known RRATs (Dencva et al. 2009; 
Keane et al. 2010; Burke-Spolaor & Bailes 2010), only 
one has a clear association with an X-ray counterpart, 
suggestive of a magnetar-like object (McLaughlin et al. 
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2007; Lyne et al. 2009). However, since RRATs are an 
observationally defined class, it may be heterogeneous. 
Many more multiwavelength associations are needed be- 
fore we can understand their underlying physics. 

Large, single-dish telescopes have pioneered the study 
of RRATs and other fast radio transients through their 
high sensitivity and relatively simple signal processing 
requirements. Unfortunately, their design also limits 
their ability to survey efficiently, since survey speed 
scales as N% D 2 , where A^ a is the number of antennas 
used, and D is their diameter (Cordes 2008; Law et al. 
2011). More fundamentally, the arcminute-scale resolu- 
tion of a large single-dish telescope is too coarse to find 
unique multiwavelength counterparts on arcsecond scales 
(Kaplan et al. 2009). Pulsar timing can obtain very pre- 
cise localizations, but requires months or years of moni- 
toring. Single-dish telescopes are also susceptible to ra- 
dio frequency interference (RFI), compromising some of 
their most exciting potential discoveries (Lorimcr ct al. 
2007; Macquart et al. 2010). 

Interferometers, being composed of a distributed set 
of antennas, can overcome these limitations. They si- 
multaneously have a large field of view and arcsecond 
resolution. Since interferometers are distributed over a 
large area, they are less susceptible to RFI and can lo- 
calize interfering sources. However, pushing interfero- 
metric approaches down to subsecond time scales has 
seemed impractical due to the massive data through- 
put and computation needed to support image-based 
analysis. Novel signal processing concepts, such as the 
"8gr8" mode at the Westerbork Synthesis Radio Tele- 
scope (Janssen ct al. 2009), "V-FASTR" (Wayth et al. 
2011), "Fly's Eye" (Siemion et al. 2012), or pulsar gat- 
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ing (Briskcn et al. 2000), have opened access to more of 
the information available to interferometers. However, 
these concepts either require a priori information or are 
not based on visibilities, so they they lose some ability 
to blindly survey and localize transients. 

Here wc describe the first practical application of very 
high time-resolution observations with an interferometer, 
recording VLA visibilities every 10 ms, and using closure 
quantities to search for radio transients. Our goal was to 
improve the localization precision of RRAT J0628+0909 
from arcminutes to arcseconds to facilitate a search for 
multiwavclcngth counterparts. Discovery observations 
with the Arecibo Observatory measured a pulse width 
of 10 ms, peak brightness of 85 mJy, and localization un- 
certainty of 3'.5 (Cordes et al. 2006; Deneva et al. 2009), 
showing that it is accessible to our new VLA observing 
mode. Our observation marks the first use of the full 
VLA at this time scale and the first interferometric de- 
tection and localization of a fast radio transient. This 
also provides the first blind demonstration of a tran- 
sient detection algorithm based on interferometric clo- 
sure quantities (Law & Bower 2012) and a second exam- 
ple of how to implement a general-purpose, fast corre- 
lator (Law ct al. 2011). We use the arcsecond localiza- 
tion of the transient to exclude its association with any 
cataloged multiwavclcngth counterpart, which excludes 
its association with a luminous, young, isolated neutron 
star. 



TABLE 1 
Observing Configuration 



2. OBSERVATIONS 



2.1. 



Developing a New Observing Mode 

Our observations required developing a nonstandard 
observing mode for the Wideband Interferometric Digital 
ARchitecture (WIDAR), the digital system at the heart 
of the VLA (Carlson & Dewdney 2000; Perley et al. 
2011). To produce visibilities, the correlator architec- 
ture uses an "XF" design, where the cross-correlation of 
antenna voltages is done in custom FPGA hardware and 
the Fourier transform stage is performed in a commod- 
ity compute cluster known as the "Correlator Backend" . 
Through a series of tests in early 2012, we tested the 
data throughput of WIDAR with a focus on correlator 
configurations that could use the full array on time scales 
much faster than 1 s. 

In normal observing, the integration time out of the 
correlator is made as long as possible to minimize data 
volumes. The length of the integration is constrained 
by the requirement that the baselines not change signif- 
icantly during an integration; that the atmosphere and 
the instrument remain stable; and that the visibilities 
can be flagged to remove RFI and other transitory ef- 
fects, without removing too much useful astronomical 
data. For the VLA, this normally means dump times of 
~ 1 second, which for a typical observing mode (cross- 
correlations for 27 antennas, two [128 MHz] subbands, 
and all four polarization products) gave total data rates 
of 1.4 MB/s. To detect individual pulses from RRATs 
like J0628+0909 requires dumping the data 100 times 
faster (10 ms), which in turn required removing several 
bottlenecks in the correlator system. The biggest speed 
improvement came from changes in the Correlator Back- 
end software, which processes the visibilities and writes 
them out to disk; in particular, changes to the implc- 



Paramctcr 
Antennas 
Configuration 
Longest baseline 
Channels per subband 
Subband bandwidth 
Central frequencies 
Polarizations 
Integration time 
Data rate 



Value 
25, 25 m diameter 
C 

12700 A rj 3000 m 
64 

128 MHz 
1312, 1440, 1568 MHz 
2 

10 ms 
90 MB s" 1 



Direction (J2000) 
Duration 



Target Eield~ 



(6h28m33s, +9°9'00") 
16 minutes 



mentation of metadata and a reduction in the thread 
count. These and other improvements allow WIDAR 
(in some constrained correlator configurations) to achieve 
data rates as high as 140 MB/s. 

2.2. J0628+0909 

On 23 April 2012, we used the new mode to observe 
J0628+0909. Observations used three subbands with 25 
antennas and two polarization products. For our 10 ms 
integration time, the total data rate was 90 MB s _1 . 
The three, 128-MHz subbands had central frequencies 
of 1312, 1440, and 1568 MHz. The highest of these was 
strongly affected by RFI and was ignored in this analysis. 

We observed J0628+0909 for a total on-source time 
of 16 minutes. We observed J0632+103 as a complex 
gain calibrator and 3C 147 as a absolute flux density 
calibrator 2 . All observations used an integration time of 
10 ms. The total data volume for the observation was 
about 200 GB. 

After the observation, we discovered that the written 
data set was not complete. A random set of visibility 
spectra (all channels for a given subband, baseline, and 
integration) were written to disk as 0. The only trend 
seemed to be that the frequency of appearance of zeros 
increased with time to a maximum of roughly 50% of all 
visibility spectra. This trend suggested that the corre- 
lator was not able to write data to disk as fast as our 
observing configuration required, despite our earlier suc- 
cess. We flagged these dropped data, which reduced our 
sensitivity by roughly 30%. 

3. ANALYSIS 
3.1. Flagging and Calibration 

The data were calibrated using the Common Astron- 
omy Software Package (CASA) version 3.3. Data identi- 
fied as bad by the VLA on-line system were flagged, in 
addition to visibilities that were exactly equal to zero. 
Two channels with persistent RFI were flagged, but no 
time-dependent RFI flagging was done, since it is difficult 
to distinguish between RFI and astronomical transients. 
Since the bispectrum transient detection algorithm can 
distinguish RFI and transients, RFI flagging is helpful 
but not required. 

As described below, our transient detection algorithm 
is based on the bispectrum, which is immune to gain 

2 Using the NRAO "Perley-Butler 2010" flux calibration stan- 
dard. 
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phase errors. Gain amplitude calibration is needed, since 
the algorithm uses the relative values of all bispectra to 
measure spatial structure of transients. For a similar 
reason, we assume that all baselines have similar sensi- 
tivity, which is generally true for VLA observations at 
these frequencies; in the future, data weighting can be 
used to remove this assumption. Bandpass calibration 
is needed in order to coherently average the data in fre- 
quency. As such, the primary goal of calibration is to 
correct for bandpass, gain amplitude, and flux scale. 

In our case, the complex gain calibrator data were not 
useful due to human (CJL) error. We instead used the 
flux calibrator, 3C 147, for all calibration. Gain solutions 
were well behaved and images of the calibrated data to- 
ward 3C 147 showed a point source with the expected 
flux density. Since 3C 147 is more than 40° away from 
our target, it is not an ideal gain calibrator. It is reason- 
able to expect differences in gain phase relative to the 
target field due to their different paths of propagation. 
While gain amplitude calibration should transfer reliably 
between 3C 147 and our target field, imaging and astrom- 
etry are not guaranteed. We discuss these issues in detail 
in §4.1. 

3.2. Transients Search 

One of the biggest challenges of this project was to ef- 
ficiently search a large data volume for transients. Our 
strategy was to use a simple, efficient algorithm to iden- 
tify candidate transients that were followed up in a more 
detailed analysis. For a first pass, we used a transient 
detection algorithm based on the bispectrum, the prod- 
uct of three visibilities from baselines that form a closed 
triangle 3 (Cornwell 1987; Kulkarni 1989; Ro gers et al. 
1995). As described in Law & Bower (2012), the bis- 
pectrum is effective for transient detection because it is 
computationally efficient, sensitive to sources over a wide 
field of view, and immune to antenna-based gain phase 
errors. The standard deviation of bispectra is also a sim- 
ple way of quantifying whether a source is point-like (as 
expected for fast transients) or extended; this functions 
as an efficient RFI rejection criterion. 

After calibration, all analysis was done in a custom 
Python radio interferometry data analysis package called 
tpipe 4 . tpipe defines pipelines for analysis and visu- 
alization of visibilities with a focus on finding tran- 
sient sources. The package uses radio interferometry li- 
braries CASA (http://casa.nrao.edu), miriad-python 
(Williams et al. 2012), and aipy (Parsons & Backer 
2009). Its design is flexible enough to read in Measure- 
ment Set or Miriad format data for either integration- 
based or dispersion-based transients searches. 

Since the dispersion measure (DM) of J0628+0909 
is roughly 88 pc cm -3 (Deneva et al. 2009), we used 
a dispersion-based pipeline for our transients search. 
Across our lowest-frequency band, this DM introduced 
a time delay of roughly 40 ms, or four times our inte- 
gration time. Failing to correct for that delay would 
have dramatically reduced our sensitivity. So we used a 
dispersion-based search that dedispersed visibilities be- 
fore forming bispectra, images, or a phased-array beam. 

3 The phase of the bispectrum is another well-known closure 
quantity: the closure phase. 

4 See http://github.com/caseyjlaw/tpipe. 



A dispersion track is defined for a given time, DM, and 
pulse width. In all analysis presented here, we assume 
a pulse width equal to the integration time of 10 ms, as 
seen in previous detections (Deneva et al. 2009). 

The algorithm searched for transients by differencing 
visibilities in time. It started by dedispcrsing visibilities 
for all 300 baselines in our 25-antenna data set. Each set 
of dedispersed visibilities then had a mean background 
subtracted in time. The background was measured for 
each baseline and channel over a sliding window of 20 ms 
on either side of the dispersion track. The total back- 
ground duration used was 40 ms, but this was spread 
over a somewhat larger window in time due to disper- 
sion. To discern spectral structure in candidate pulses, 
each search was repeated for assumed DMs of 44, 88, and 
132 pc cm" 3 . 

After forming the background-subtracted, dedispersed 
visibilities, a bispectrum was calculated for all closed 
triangles in the array. For n a antennas, there are 
n a {n a — 1) (n a — 2)/6 bispectra, all of which contain in- 
dependent information (Kulkarni 1989). For our obser- 
vation, we calculated 2300 bispectra for each integration, 
DM, and polarization. After forming the bispectra, we 
took the mean over both polarizations to effectively pro- 
duce a Stokes I bispectrum per integration and DM. 

Because the visibilities were averaged in frequency 
before forming the bispectrum, beam smearing (chro- 
matic aberration) decreases the sensitivity to sources far 
from the field center. The FWHM of the resulting de- 
lay beam is ^dciay ~ 4(i//BW)0 s b, where v is the ob- 
serving frequency, BW is the bandwidth used in form- 
ing the bispectrum, and 8 s b is the synthesized beam 
size (Thompson et al. 2001; Taylor et al. 1999). Our 
search used a bandwidth of 256 MHz and a beam size 
of roughly 15", producing a delay beam of 5'. This is 
large enough to cover the localization uncertainty of 3'.5 
for J0628+0909 (Deneva et al. 2009), so the algorithm 
could be applied as described above 5 . 

We searched for astronomical transients using both the 
real part of the mean and the standard deviation of the 
bispectra at each time and DM. The real part of the 
mean bispectrum detects any anomalous event, since it 
is proportional to the cube of the brightness of an event 
(both transients and RFI). The standard deviation over 
all bispectra in the array grows roughly as the square of 
the event brightness, if it is a point source. RFI has a 
different effect on the standard deviation of the bispec- 
tra, since it tends to be near the telescope or affect each 
baseline differently. The net effect is that RFI looks like 
a spatially extended transient, which has a large stan- 
dard deviation over bispectra (Law & Bower 2012). So 
to select transients consistent with a celestial source, we 
required (1) a 5<r deviation of the mean bispectrum and 
(2) a standard deviation equal or less than that expected 
for a point source. For this and other datasets, this ap- 

5 In cases where the delay beam is too small, one could form 
multiple delay beams or form the bispectrum from an average of 
frequency segments that have a larger delay beam. In the latter 
case, for JV frequency segments, the noise per segment scales as 
TV 1 / 2 , but the SNR for the bispectrum scales as the SNR cubed, 
or N~3/ 2 (at low SNR; Rogers et al. 1995). Averaging bispectra 
over segments benefits from another factor of N 1 ' 2 so overall the 
bispectrum loses sensitivity linearly with the number of segments 
used. 
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Fig. 1. — (Top:) Bispectrum light curve for roughly 4 seconds of 
1 subband of VLA data toward J0628+0909. The red star shows 
the pulse detected by the algorithm. (Bottom:) Mean versus the 
standard deviation of bispectra across the array for each integra- 
tion. The standard deviation is normalized by an estimate of the 
noise per baseline cubed. The solid line shows the theoretical pre- 
diction for a point source with a range of brightnesses; a single, 
transient point source will appear on that line. The dashed line 
is the same as the solid line, but shifted to the right by 20%. To 
tolerate some variance due to noise, we consider any point left of 
the dashed line to be a candidate. Both panels are scaled to show 
the transient, so RFI lies far above (in the top panel) or to the 
right (in the bottom panel). 

proach produced a list of candidate events that can be 
easily followed up manually. 

Our first step in following up any candidate event was 
to look for associated RFI. Some RFI may have point- 
like spatial structure (e.g., from satellites), but such RFI 
almost always has narrow (< 5 MHz) spectral structure. 
We searched for RFI by creating a spectrogram for times 
near the candidate. We did this by summing cross cor- 
relations to effectively form a phased-array beam at the 
phase center, assuming a phase-calibrated array. The 
vast majority of RFI was evident in a spectrogram like 
this. A second test of any candidate event was to image 
the dispersion trial and look for a point source. Finally, 
we moved the phase center to the location of any can- 
didate source seen in an image, formed a spectrogram, 
and looked for the characteristic quadratic sweep of a 
dispersed pulse. 

4. RESULTS 

4.1. Event Detected with Bispectrum 

We used the bispectrum transient detection algorithm 
to search 16 minutes of VLA data with two, 128-MHz 
subbands toward J06284-0909. With 10 ms integrations 
and three trial DMs, we searched roughly 2.8 x 10 6 dis- 
persion trials and found 14 candidate events. Of those 
candidates, all but one were traced to RFI or otherwise 
corrupted data. Figure 1 summarizes the bispectrum 
algorithm during the candidate event when using only 
the lowest subband. It is detected by the bispectrum 
with roughly 16a significance 6 and has a standard devia- 

6 For convenience, we quote the apparent SNR measured in the 
mean bispectrum. In this case, the mean bispectrum is wider than 
a Gaussian, so the significance is overstated by a few percent. 



Fig. 2. — Image of the dispersion trial with the candidate tran- 
sient detected by the bispectrum algorithm and shown in Figure 
1. The image has a beam size of 21"xl5"; the colorbar shows the 
flux scale in units of Janskys. A point source of significance 13<x is 
seen within 1' of the previous estimate for the location of RRAT 
J0628+0909. 

tion consistent with the theoretical prediction for a point 
source. 

For the candidate transient not clearly associated with 
RFI, we ran a series of tests to confirm the signal and 
show that it was associated with J0628+0909. First, the 
bispectrum search was repeated for data from the first 
and second subband independently. Similarly for polar- 
ization, we ran the algorithm for the RR and LL cor- 
related products independently. The transient candidate 
was detected in all of these cases, showing that it is spec- 
trally broad and not strongly circularly polarized. 

Next, we imaged the background-subtracted dispersion 
trial associated with the candidate to confirm that it is 
point-like, as indicated by the standard deviation of the 
bispectra. Figure 2 shows an image of the first subband 
of the VLA data for this candidate. The image contains 
a single point source with brightness of 215 ± 17 mJy 
beam -1 , or 13c significance. 

For an observation with 25 antennas, an event detected 
with apparent 16a significance in the bispectrum should 
have roughly 15a significance in the image plane, assum- 
ing good phase calibration (Law & Bower 2012). The 
slight reduction in the significance of detection in the 
image plane is consistent with the fact that the image 
noise is 20% higher than expected (considering dropped 
data) . These two points suggest that the gain calibration 
was good to roughly 20%, consistent with the reasonable 
image quality. We note that imaging with two subbands 
shows imaging artifacts indicative of stronger phase er- 
rors, so we focus our subsequent analysis on the lowest 
subband, covering frequencies from 1248 to 1376 MHz. 

We fit an elliptical Gaussian representative of the syn- 
thesized beam to the point source to measure its loca- 
tion and uncertainty. The point source is offset from 
the J2000 phase center of (6h28m33s, +9°09'00") by 
(48'.'1, 14'.'9), which places it at (6d28m36.25s, 9°9'14'/9). 
This offset is within the localization uncertainty of pre- 
vious detections (3'.5), which clearly associates it with 
J0628+0909. It is also within the delay beam of our 
search (5'), which shows that it is the same pulse seen 
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Fig. 3. — Spectrogram showing the dispersed pulse from RRAT 
J0628+0909 shown in Figure 1. The spectrogram is generated for 
a phased-array beam pointed at the transient seen in the image in 
Figure 2. The color scale is in units of Janskys. 

by the bispectrum algorithm. The uncertainty in the 
Gaussian centroid in (RA, Dec) is (l'/3, l'/5). 

To test the astromctric accuracy of the data, we mea- 
sured location of a steady source elsewhere in the field. 
We imaged NVSSJ062859+091149 in a series of 10 1- 
s snapshot images around the time of the bright pulse 
from J0628+0909. The cataloged source location is 
(6:28:59.99, +9:11:49.7) known to a precision of roughly 
1". The mean position measured in our snapshot images 
is (6:28:59.93, +9:11:49.5) with a standard deviation of 
(l'.'O, 0'.'6) and confirms its astrometric accuracy. Adding 
this uncertainty in quadrature to the source centroid er- 
ror, we find an absolute centroid precision of (l'/7, l'/6). 

As final confirmation of the detection, we formed a 
beam on the measured location of the pulse and gener- 
ated a spectrogram. Figure 3 shows the spectrogram of 
the first subband toward the transient. The transient is 
clearly a broad-band pulse with an unresolved width of 
10 ms. The pulse spans five integrations with the char- 
acteristic quadratic sweep expected of dispersion. 

Figure 4 shows the apparent pulse brightness in a 
phased-array beam as a function of assumed DM and 
pulse time. This kind of analysis mimics traditional, 
single-dish pulse searches, which rely on the dispersive 
delay to distinguish an astronomical source from RFI. 
We find a peak pulse brightness for DM near 90 pc cm , 
consistent with previous detections. 

4.2. Events Detected with Phased Array 

With a priori constraints on the location of 
J0628+0909, we searched for pulses with a phased-array, 
dispersion-based search. We expected to detect more 
pulses with a phased-array beam than with the bispec- 
trum, since the bispectrum sensitivity falls off as the cube 
of source brightness for SNR per baseline less than 1 
(Rogers et al. 1995). For our single VLA subband data 
with 10 ms integrations, this happens for a brightness of 
roughly 280 mJy beam -1 . 

Since the data are largely unflagged, a phased-array 
beam is highly susceptible to RFI. We filtered candidate 
transients by finding events with peak flux for an as- 



FlG. 4. — Pulse SNR seen in phased-array data for a range of dis- 
persion trials toward RRAT J0628+0909. Circles mark dispersion 
trials with SNR greater than 5 and scale linearly with SNR. 

sumed DM near 88 pc cm" 3 . As with candidates from 
the bispectrum search, we rejected candidates associated 
with RFI seen in spectrograms. This search revealed two 
further candidate pulses with peak phased-array beam 
SNR of 6 and 7. Imaging dispersion trials at their peak 
SNR shows point sources at the same location as the 
pulse detected with the bispectrum algorithm. All told, 
our VLA observations find three pulses with brightnesses 
of 75, 86, 215 mJy beam" 1 in 10 ms. These flux densities 
may be supressed by up to 20% due to phase calibration 
errors. 

With three pulse detections, we can use the relative 
pulse arrival times to estimate the underlying rotation 
period of the RRAT. Assuming that each pulse arrives 
from a random rotation of the neutron star, the period 
is at most the greatest common divisor of their relative 
arrival times (McLaughlin ct al. 2006). We find that a 
period of 1.24148 s is most consistent with the rotation 
period of the neutron star, confirming earlier work. 

4.3. Using the Improved Localization 

Our arcsecond localization of J0628+0909 makes it eas- 
ier to search for counterparts. No compelling counter- 
part was found in a search with the Simbad, NED, and 
HEASARC systems. A search of X-ray data archives did 
not find any Chandra observation that could potentially 
detect thermal emission from the neutron star. 

A search of the Sloan Digital Sky Survey (SDSS-III, 
data release 9; York et al. 2000) found one potential as- 
sociation for the RRAT. Figure 5 shows a SDSS-III i'- 
band image with the large ellipse showing the la un- 
certainty on our best-fit location for J0628+0909. A 
red stellar source, SDSSJ062836.22+090917.5, is offset 
by 2cr from the radio source centroid. We use the offset 
of these two sources and the density of sources brighter 
than SDSSJ062836. 22+090917.5 to estimate the proba- 
bility of random association (e.g., Croft et al. 2011). For 
a i'-band source density of 4.4xl0 -3 and a radio source 
offset of 3'.'5, we expect a 15% chance of false association, 
suggesting that the association is not significant. 

Independent of our VLA observation, the team that 
discovered J0628+0909 has been monitoring it (Nice et 
al 2012, in prep). In Arecibo observations spanning over 
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Fig. 5.— SDSS i'-band image of the region around J0628+0909. 
The large ellipse marks the localization uncertainty for the bright- 
est radio pulse we detect from J0628+0909. This includes uncer- 
tainty in the absolute astrometry described in §4.1. The plus shows 
the location of SDSSJ062836. 22+090917.5, the nearest SDSS coun- 
terpart. The small ellipse shows a new localization of J0628+0909 
from an observing campaign with the Arecibo Observatory (Nice 
et al 2012, in prep). 

three years, they have collected thousands of pulses that 
allow a precise localization of the RRAT using tradi- 
tional pulsar timing techniques (Hobbs ct al. 2006). Us- 
ing these observations, they have localized J0628+0909 
to (6:28:36.1832±0.0025, +9:9:13.93±0.18). As shown in 
Figure 5, the timing solution is within the la uncertainty 
in our localization, confirming the image-based analysis. 

By excluding an SDSS counterpart to J0628+0909, we 
can put an upper limit on its optical luminosity. The 
typical limiting magnitude is i'= 21.3 mag, which corre- 
sponds to an apparent flux density limit of 1.1 x 10~ 5 Jy 
(Fukugita et al. 1996). The dispersion of pulses toward 
J0628+0909 is 88 pc cm -3 , implying a distance of 2.5 
kpc (typically better than 30% accurate; Cordes & Lazio 
2002). At that distance, we expect extinctions of Ay = 
0.98 mag and A v = 0.62 mag (Drimmcl ct al. 2003; 
Cardelli et al. 1989). Together, these effects imply a up- 
per limit on the i' luminosity of 1.1 x 10 31 erg s _1 . It 
is difficult to compare this to the optical/infrared prop- 
erties of isolated neutron stars, since observations - 
and the spectral and temporal behavior of neutron stars 
themselves — are heterogeneous (Rca et al. 2010). How- 
ever, the SDSS i' limit does exclude luminous, young neu- 
tron stars (e.g., the Crab; Eikenberry et al. 1997), which 
have flat infrared spectra and can be orders of magnitude 
brighter than the SDSS limit. 
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5. CONCLUSION 

We have demonstrated a new observing mode for the 
VLA that can write 10 millisecond visibilities, opening 
access to this exciting timcscalc to image and visibility 
analysis techniques. We used this new observing mode 
to find and image individual pulses from the poorly- 
localized RRAT J0628+0909. This is the first blind 
demonstration of our pulse detection algorithm based on 
the interfcromctric closure quantity known as the bispec- 
trum and the first use of visibilities to detect and localize 
a fast radio transient. 

We measure a new position for J0628+0909 of 
(6d28m36.25s, 9°9'14'/9) (J2000) with an accuracy of 
roughly l'/6. Our interferometric localization is two or- 
ders of magnitude better than the localization made in 
its discovery observation. It is also consistent with new, 
sub-arcsecond localization made via a traditional tim- 
ing analysis with Arecibo data spanning more than three 
years. Our localization for J0628+0909 excludes its asso- 
ciation with any cataloged multiwavclcngth counterpart. 
SDSS imaging places an upper limit on the RRAT op- 
tical luminosity that excludes association with a young, 
luminous, isolated neutron star. 

This paper presents the first real-world application of 
new signal processing and algorithms for using radio in- 
terferometers to the study millisecond transients. More 
than just pushing traditional slow transient techniques 
to fast time scales, our study of J0628+0909 shows how 
fast interferometric imaging encompasses and extends 
single-dish techniques. Future iterations can expand this 
concept by implementing real-time, all-time processing, 
opening the latent survey ability of radio interferometers. 
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